The Shubnikov-de Haas effect in multiband quasi-two dimensional 

metals 



I. O. Thomas^, V. V. Kabanov^ and A. S. Alexandrov^ 
^Department of Physics, Loughborough University, 
Loughborough, United Kingdom, LEll 3TU 
Josef Stefan Institute 1001, Ljubljana, Slovenia 

Abstract 

We analyze the behavior of the longitudinal conductivity cTzz in a field perpendicular to the 
highly conducting plane of a quasi-2D multiband metal in the case of closed system and open 
systems where r is fixed, r oscillates due to intra-band scattering and r oscillates due to inter- 
band scattering. In all but one case, we find that there are also mixing frequencies present - 
however they exhibit different qualitative behaviors, as befits their different origins, and in the 
case of inter-band scattering in an opensystem, may in fact be absent in the dHvA oscillations of 
that system. 

PACS numbers: 72.15.Gd, 75.75.+a, 73.63.-b 



1 



I. INTRODUCTION 



The use of magnetic quantum oscillations [l| as a means of analyzing the Fermi surfaces 
of low- dimensional materials has received interest in recent years, particularly with respect 
to organic compounds such as the metallic varieties of the charge transfer salts BEDT-TFF 
(see reviews by Singleton 2|] and Kartsovnik [3]). One interesting feature of such systems 
is that the magnetic quantum oscillations can render measurable an apparent difference 
between the canonical (closed) and grand canonical (open) ensembles (in which the density 
of electrons Uf. or the chemical potential are fixed, respectively); in the former case the 
chemical potential /i will oscillate, and this will have a noticeable effect on the behavior of 
the magnetic oscillations of a two- or quasi-two-dimensional system. 



In two dimensional multi-band metals, frequency mixing in the de Haas-van Alphen 
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potential magneto-oscillations has been predicted 14 



effect resulting from chemica 
analytically and numerically 



B,y,0,y,y,y 



11| and observed experimentally 



studied 
in 



^Qii 



closed systems. The theoretical description has been generalized to the more realistic quasi- 
two-dimensional (quasi-2D) case [ij]. Recently, it was proposed that finite amplitudes due to 
mixing between extremal cross-sections of the warped Fermi surface of a single band closed 
<,ua.-2D ^ysteo, ate be observed fl; unlike the — exa.,ple.. th. particular 

mixing effect cannot be replicated by magnetic breakdown P, |16[ . 

Attention has recently turned to the Shubnikov-De Haas (SdH) oscillations of the longi- 
tudinal conductivity of quasi-2D single-band metals, with quantum transport theory being 
used by Champel and Mineev 17(1 to examine the ultra-high field, 2D limit, and by Grigoriev 



[isl to examine the intermediate (quasi-2D~ 
physics of the former case remains obscure 



field limit. While much regarding the detailed 



191], analysis of the latter case while taking into 



account scattering in the self-consistent Born approximation with t ^ T^, where is 
the Dingle temperature and t the interplane transfer integral, has resulted in a proposed 
explanation of a few oddities in the observed oscillations - namely, the existence of slow 
oscillations and the field dependent phase shifts of the beats of the conductivity. 

In what follows, we examine the behavior of the longitudinal conductivity in a field 
perpendicular to the highly conducting plane of a quasi-2D multihand metal. We begin 
by deriving an expression for the conductivity of a two-band metal in the abscence of any 
scattering or chemical potential oscillations in Section III A[ We then (Section IIIBI) derive 
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an expression for the self-energy of a system in the Born approximation where there ex- 
ists scattering between two bands, which we will have need of in some of our following 
calculations. 

In the section following that, we calculate the effects of chemical potential oscillations 
and different forms of scattering on the system. Firstly, in Section IIII Al we consider a 
closed system with a fixed relaxation time r, then in Section IIIIBI an open system with 
oscillations in r arising from intra-band scattering (this is merely a two-band generalisation 



of Grigoriev's 18| result), in Section IIII CI a closed system with oscillations in r arising 
from intra-band scattering, in Section IIIIDI an open system with oscillations in r arising 
from inter-band scattering, and finally in Section IIII El a closed system with oscillations 
in T arising from inter-band scattering . In all cases other than the second we discover 
some measure of frequency mixing, though the precise behaviour of such mixing varies 
according to its origin. Interestingly, as argued at the beginning of Section IIIIl one would 
probably not expect frequency mixing due to inter-band scattering to be very significant in 
the dHvA effect; that it is present in open systems in the SdH effect seems to indicate that 
one distinction between open and closed low-dimensional thermodynamic systems is slightly 
blurred in the case of conductivity. 

II. SDH EFFECT WITH FIXED RELAXATION RATE 
A. General expression for longitudinal conductivity 

A quasi-2D metal in a longitudinal magnetic field has the following energy spectrum 
,3: 

ea{n, kz) = uJain + 1/2) — 2t cos^k^d), (1) 

where a labels the band, d is the distance between layers, n labels the Landau level, k^ is 
the momentum in the perpendicular direction and cJq, = eB/rria is the cyclotron frequency, 
ma being the effective mass for that band, and h = c = kB = 1- 



We may calculate the conductivity from the Kubo formula 201] 

t2 /"OO 



a.. = — / rfe[-<(0]Tr m - H)v-M ' H)vz], (2) 

with Vz being the velocity operator for the 2;-direction, n'p being the derivative of the Fermi 
function, the trace being taken over all the single particle states /5 = {n, kz, ky, a} = {7, a} 



and the spin, V is the volume of the system, and the bar denotes the averaging over the 
random distribution of impurities in the sample. We may, if the impurities are point-like, 
their concentration is low and t ^ To, neglect the vertex corrections to the average and 



replace 6{C, — ep) with ^G^{ej3,^) /-n, where 

is the retarded Green's function and S^(^) = Lp{^) — iAp^C,) is the retarded self-energy 
(L^(^) and A/3(,^) being real), to obtain 

a 7 

This is simply the initial conductivity formula of references [itI. [isl. only now generalized 
so that it sums over multiple bands. Restricting our interest to the case of two bands with 
different masses mi and m2, and since t ^ T^, we may use the approximation S^(^) ^ S(^), 
makin g th e calculation analytically tractable). Then we may generalize the single-band 
result 171, [l8| as 



Here, J\{x) is the first-order Bessel function, ^* = ,^ — L(^), L(,^) being the oscillating 
portion of -L(^), = eB j^dn and 



One should use the expansion Ji{kx)/k = x/2 for the /c = harmonic. 

To begin with, let us assume that the quantum oscillations in the self-energy are 
suppressed, with the result that we can ignore the oscillatory contribution to ^* and 
A(^) = (2r)^^ = ttTd, where r is the mean scattering time and To is the associated Dingle 
temperature. This is justified in the presence of a field and size independent reservoir of 
states l4, llTI], for example, which suppresses both chemical potential oscillations and any 



oscillations in the self-energy. Having made this assumption, we integrate over ^ (using the 
delta- function-like behavior of n'p{^) near the Fermi energy at small T to obtain the first 
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term, and cos(ax) cosh '^{x)dx = aw/ sinh(a7r/2) to obtain the damping factor for the 
oscillatory terms) and so acquire: 



a 



k=l 



k 



' ^ fA7rkt\ fuJa , 27rkTD\ ( 



27TkfI 







2n^kT 
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where a = {e^dt^M) / {h'^Td), M = mi +m2, i?„ = exp{—2TT'^Tij/uJa) is the Dingle reduction 
factor, and Q{y) = y/ sinh(?/) is the usual Lifshitz-Kosevitch reduction factor. If Ant > uj^ 
we can use the asymptotic 

/ A'Kkt\ 
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and so obtain 
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where 



A1 



-Rn 
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We can see then that the conductivity oscillates with two fundamental periods, one 
corresponding to the first band and one to the second. In the quasi-3D limit one will 
observe the splitting of each fundamental in the Fourier Transform of the conductivity into 
a pair of peaks (a low- frequency 'neck' peak and a high-frequency 'belly' peak), which is 
evidence of the slight warping of the Fermi-surface such that there are two extremal orbits 
present [l|. 



B. Self-energy with interband scattering 

If scattering between bands is possible, then the self energy of the particles in any given 
band will contain contributions from all other bands in addition to that due to intra-band 
scattering. In that case, in order to obtain the conductivity in the self-consistent Born 
approximation in the two band case, we must include the impurity diagrams from both 



bands in S/3(.C) |2l| : 



E^(0 = W5^5^G^(6„,„0, (11) 



a 7 
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where W is the square of the scattering amphtude, which is proportional to the impurity 
concentration and set equal to the same constant in every band for simplicity. It also includes 
a factor of 2 from the summation over the spin. Following reference [l8|] we find that 



Vnia 



+ ( 1 + 2|;(-1)'J„ (^^) COS fH^^ ) R^kX) 



(12) 



From this, we can see that there is an oscillatory contribution from the real part of the 
self-energy that must be taken into account in our calculations, and a slowly varying part 
A{^) which may be ignored. Using our symmetric approximation (equation ffTTl) ) and setting 
WVM/2TTd = ttTd (as the average Dingle temperature is simply related to the average value 
of the imaginary part of the self-energy due to scattering) we may write the imaginary and 
oscillating real parts of the self energy as follows: 



A(0 = tcTd 



1 + 2M-^ Yl E(-l)'^o (^) cos ) R^{k, 

a k=l \ a / \ a 



m = 27rr^M-iXIm„f;(-l)Vo sin ) R^{k,0- 



(13) 



(14) 



Equation f|T3|) is a non-linear equation for A(,^) which can be solved approximately in the 
strong damping limit Ra{k,^) <^ 1, as we shall see in section UlI D[ 



III. THE EFFECTS OF SCATTERING AND OSCILLATING ON THE CON- 
DUCTIVITY 



In this section, we discuss the effects of chemical potential oscillations and of various 
kinds of scattering on the conductivity. It should be noted that this analysis takes place in 
a region where there is strong damping of the amplitudes by the factor Ra - we do not work 
in the region where the processes giving rise to mixing frequencies are at their strongest, 
and that given that the effects we are interested in are of the second order in R^, they may 
be quite small. However, it has been noted (by comparison to numerical results) in the case 



of dHvA oscillations 



15| that a reasonable level of accuracy may be maintained even if 



one allows the Dingle damping factor to tend towards unity, assuming that one is in the 
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appropriate limit regarding the behaviour of /i. This might also be the case with respect to 
the SdH oscillations where mixing due to inter-band scattering is present. 

Indeed, it would be suprising if our analysis of those systems were not even qualitiatively 
correct (with regards to the presence of frequency mixing at the very least) outside of the 
region where ^ 1, since the amplitudes contain a denominator A(^) which becomes 
more important as it becomes smaller and Ra therefore approaches unity. Interestingly, this 
would not be the case with regards to dHvA oscillations - as can be seen from equation (68) 
of 



22l |. A(^) enters only through the Dingle damping factor and the cosine function. This 
entails that, as Ra approaches 1, the oscillating portions of A(^) become negligible - it is 
doubtful, therefore, that one would be able to observe frequency mixing due to scattering 
in the dHvA oscillations of any but the dirtiest systems. 

In cases where slow oscillations of the conductivity exist, it should be noted that it has 
been observed 



23( 1 that macroscopic inhomogeneities in the sample increase the damping of 
the fast oscillations over and above that expected by simple scattering; this can be modeled 
by replacing the Dingle temperature in the Ra factors of those oscillations with a new, 
larger Dingle temperature that will enhance the damping of their amplitudes. In order 
to simplify our discussion, though, we ignore this complication for now, though it should be 
taken into account in any physical measurement. 

A. Canonical Ensemble: oscillating fixed r 

In closed systems, the particle density is fixed and we work in the Canonical Ensemble. 
As a result of this, the chemical potential /i of the system may oscillate. In three dimensions 
these oscillations are negligible, however, as mentioned previously, in the case of multiband 
2D and quasi-2D metals with two bands one discovers additional harmonics corresponding to 
frequency mixing between bands. In this section, we discuss the effects of these oscillations 
on the conductivity of a system in which scattering effects are suppressed. This situation is 
physical if the oscillations in r are suppressed due to a large amphtude of scattering from 
quantised bands to the resevoir, which is consistent with the existence of a small enough 



resevoir density of states 
Following reference 



or the oscillations in yU to remain significant, 
we divide the chemical potential into an oscillating part fi and a 
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non-oscillating part /iq: 

H = fio + fl, (15) 

where 

with Q being the oscillatory portion of the thermodynamic potential Q, p = pbg + "^^Pa, 
Pbg being the unquantized background density of states which we treat as being negligible, 
and Pa being the quantized density of states for the band a. Using formula (14) of reference 
jl^ and the 2D density of states pa = "^o/vr, we acquire 

,.-j:±Bts^'-^^)B(?^), (17) 



a k=l 



where 



2eB{Raf f4:7ikt\ 



= ^^Jf 



nMk 

which follows from the use of the asymptotic 



Ua f Ankt tt\ 



Let us assume that the system is strongly damped: Ra ^ 1, and we keep only the 
oscillating terms of second order in Ra or less. It follows that fi/uja <^ 1, so we may expand 
([9]) in /i and approximate po with p. We then insert (fTTIl . again approximating the argument 
of the sine function as /i ~ /io, and so obtain 

where a°l'^ represents the unmixed portion of the conductivity and cr™^ is the mixed portion. 
(T°l'^ is given by: 



a ^ V \ UJa J \UJa 



+ 



maA'^e ( 1 + CaaQ 



^a \ UJ, 



COS 



iirp 



-CaaQi^) 1, (21) 



which shows that the oscillations in ft, modify the second harmonics and also create additional 
slow oscillations due to the effect of the warping of the Fermi surface on the behaviour of the 
23], as predicted for the single band case by Grigoriev This slow oscillation 



scattermg 



is not seen in dHvA oscillations (such as those of reference 15|], for example); in those cases, 
one sees only the fast oscillations whose frequencies are determined by t and fi together, 
with no slow oscillations whose frequencies are determined by t alone. 
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FIG. 1: A graph showing the Fourier transform of the SdH oscillations from .9 < Bq/B < .95 
where the frequencies are given in terms of B/Bq, alongside plots of how the damping factors 
evolve as their respective temperatures are increased. Due to the small size of the window, the 
slow oscillations are poorly resolved. The y axis has been rescaled by a factor of 100. The following 
parameters are used: fi = 500t, {uji/ijj2) = (m2/mi) = .9, ksTo = .026t and ksT = .GOOOSt, 
setting our unit of measurement to be Bq = mit/2he, which is around 41 Tesla if t = O.OleV and 
mi = rrie. The legend in the Fourier Transform plots should be interpreted as follows: A - closed 
system, fixed r; B - open system, inter-band scattering; C - closed system, intra-band scattering; 
D - closed system, interband scattering; E - open system, intraband scattering. The legend in 
the Temperature Damping plots should be interpreted as: a - value of 0(27r^T/a;a); b - value of 
@{27t'^T /uji)@{2tt'^T /uj2) . The legend in the Dingle Damping plots should be interpreted as: I - 
value of Ra] 11 - value of R1R2', HI - value of i?^. 
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The mixed component of the conductivity is given by 



a 



M-\Ci2 + C2i)Q 



2fT 







2fT 

U2 



COS 



2-71 fl 



UJ4 



COS 



where 



1 _ 1 ^ 1 

UJ± UJ2 UJl ' 



(22) 



(23) 



and it is this term which gives rise to the frequency mixing induced by the chemical potential 
oscillations. 

The amplitudes induced by the chemical potential oscillation are given by: 



r< — 

in the case where a = a', and 
_ 2eBmaai^ 



MttH 





'Ant 


TT 






4_ 



MttH 

in the mixed (a 7^ a') case, where 



niryi 



1/2 



COS 



47rt\ 

UJ+J 



± sin 



rvrt t 



(24) 



(25) 



(26) 



and we take the positive sign in front of the sine when a' = 1 and the negative when a' = 2. 

One important feature of the mixing amplitudes in this case is that they are identical for 
both the additive and the subtractive mixing frequencies. This does not generally hold, as 
we shall see in what follows. 



B. Grand Canonical Ensemble: Intra-band scattering 



Working in the self-consistent Born approximation, let us assume that there is no scat- 
tering between the bands, and that we work in an open system described by the Grand 
Canonical Ensemble. In this case the only diagrams contributing to the self energy of an 
electron in band a will be those corresponding to intra-band scattering, and there are no 
chemical potential oscillations that could also result in a mixing of oscillation frequencies. 
This is the situation described by Grigoriev ISj. Here we simply generalize it to the case of 
multiple bands. Keeping all terms up to and including 0(i?^), we obtain: 



1 + j: ™„ ( are ( COB (?^) + f Vo= I + fl" 



(jJn 



(27) 
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where 



= 2 



77%: (1 + {ai,ay)Ra COS ( — - Y + 



(28) 



D2 = — ^Vl + (ai,„)2(i?„)2cos 7 cos 7 + 01,a 



2 I UJa ,^ , sr.. I Svrt TT 



+ (i?„)^/^ (1 + (a2,„)2) COS ( — - - 



_ CJa(-Ra)'^ 



A-Rt TT 

7 + '^50 



= arctan(ar^fj), 05^^ = arctan 



(29) 
(30) 



and Or, a is given by equation (l26l) . 



C. Canonical Ensemble: oscillating /x, intra-band scattering 

When the system is closed, we must take into account the effects of the oscillations in ix 
as well as that of scattering. We can do this by applying the proceedure outlined in Section 
nil Al to equation (1271) . assuming that -Rq ^ 1. 

Having performed the expansion in terms of /i, we may write our result as: 



Ozz = (^zz + ^zz ) 



(31) 



where 



and 
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The amplitudes due to the chemical potential oscillations are given by 



UJ4 



cos((/)i.Q,) + COS 2 



4'7rt vr 
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(32) 



• (33) 



(34) 
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FIG. 2: A graph showing the Fourier transform of the SdH oscillations from .9 < Bq/B < .95 with 
{uJi/uj2) = (m.2/mi) = .75, alongside plots of how the damping factors evolve as their respective 
temperatures are increased. Due to the small size of the window, the slow oscillations are poorly 
resolved. The y axis has been rescaled by a factor of 100. Except as otherwise noted, all parameters 
and legends are identical to those in Figured) 



when a = a' and by the following when a ^ a': 

1/2 



M-kH 



TTlr 



nin' 



A-nt 
cos I ± 



- Sm h 01, a 



(35) 



where the argument of the cosine function contains a '+' if a = 2 and a ' — ' if a = 1. 

In this case, we find that we have two terms contributing to the slow oscillations, one 
of which is temperature dependent, and that (as one might expect from section IIII Ap we 
have additional frequencies due to the mixing of bands by the oscillations in the chemical 
potential. 



D. Grand Canonical Ensemble: inter-band scattering 

In order to proceed with our calculation of the effects of inter-band scattering on the 
behaviour of the system in the Grand Canonical Ensemble, we expand ( fT3l) in the strong 
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damping limit Ra{k,^) ^ 1, where we make the approximation 



1 H — 2_,'^a'Jo y~j I — 1 
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FIG. 3: A graph showing the Fourier transform of the SdH oscillations from -9 < Bq/B < .95 with 
{u!i/ui2) = {m2/mi) = .6, alongside plots of how the damping factors evolve as their respective 
temperatures are increased. Due to the small size of the window, the slow oscillations are poorly 
resolved. The y axis has been rescaled by a factor of 100. Except as otherwise noted, all parameters 
and legends are identical to those in Figured) 

We can then expand out the cosine terms for small L{^), and find that for k = 1: 

f A'nt\ f A-n^ 



cos I 1 i?Q,(l,4) ^ cos ( ) i?a + 



[Ra) maJo COS 



H 7-rRaRamaJo COS 



(37) 



where if a = 1, 2 then a = 2, 1, and that for k = 2: 



cos I — - I i?a(2,0 ^ COS ( — 1 Ri. 

UJry / \ UJn 



(38) 



The conductivity in this approximation is given by (after the substitution of f|T3l) into ([5])): 

/47rt\ f2nC* 
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(39) 



where cTq = {e^dt^nia) / {•n'^To)- 

Substituting fl36|) . fl371) and (138|) into (l39l) . integrating over ^, replacing the Bessel functions 
with their asymptotics and gathering all the terms together, we finally obtain: 
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Here, the unmixed amplitudes are: 
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As in the case of oscillating /i, we observe mixed frequencies as well: 
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where 
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FIG. 4: A graph showing the Fourier transform of the SdH oscillations from -9 < Bq/B < .95 with 
{uji/lj02) = {1712/ nil) = -4, alongside plots of how the damping factors evolve as their respective 
temperatures are increased. Due to the small size of the window, the slow oscillations are poorly 
resolved. The y axis has been rescaled by a factor of 100. Except as otherwise noted, all parameters 
and legends are identical to those in Figure [H 



As before, the values of a^^a and 4>r,a and (psa are given by ( l26|) and ( l30l) respectively. 



E. Canonical Ensemble: oscillating fi, inter-band scattering 

Let us now treat the above system as though it were closed, and so allow fi to oscillate. 
The only terms of interest to us will come from the expansion of the terms proportional to 
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. From examination of (l28i) and (HTj) we can see that 

^ M ^' 

and it follows that the amplitudes due to the chemical potential oscillations are: 
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and that the expression for the conductivity is: 
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(50) 



In this case, all the second order terms are modified by the chemical potential oscillations. 
IV. DISCUSSION 

Figures [1] to m display for purposes of illustration the Fourier transforms of the oscillatory 
components azz of (cr^z— o")/cr given by the expressions in the previous sections. The inclusion 
of plots of the behaviour of the Dingle and temperature reduction factors in the figures 
should facilitate the extrapolation the results given here to regimes where -C 1, where 
our analytic calculations are more valid; however, it would be suprising if the qualitative 
elements of our results were not preserved even at values of Ra near to unity. For the 



purposes of reference to a real material. Cole et a/.'s [2J] measurements of the effective 
mass in the two-band system of GaAs-(Ga,Al)As hetero junctions indicate a ratio of masses 
somewhere in the region .3 ^ (1712/1711) ^ .45. 

Note that in order to clean up the signal and remove spurious 'ringing' due to the finite 
window size, the data was passed through a Hanning window prior to the operation of the 
numerical Fourier transform (see, for example, reference 25l|). 

Details of the mixing frequencies are shown in Figures [5] and [61 Perhaps the most obvious 
and interesting feature is that the signature of the two kinds of mixing is not the same: the 
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FIG. 5: Figures showing the details of the I/2 — fi\ amphtudes of (771-2/7711) values of a) 0.9, b) 0.75, 
c) 0.6 and d) 0.4 . 

amplitudes of the I/2 + /1I and the I/2 — /i| mixing are identical in the cases where the mixing 
is due to only oscillations in fi, but in cases where mixing through inter-band scattering is 
present, the I/2 — /i| mixing amplitude is less than that of the I/2 + /1I mixing (the exception 
is case C at a ratio of .9, most likely due to the small size of its amplitude and its being 
interfered with by neighboring peaks). In general, the amplitude of the mixing peaks grows 
more pronounced as the two bands become less similar (that is, as m2 becomes smaller than 
mi), at least for these values of 1712, apart from case A at a ratio of .4, where a splitting of 
the peaks reduces their size. It should also be noted from the structure of the LK factors 
in the mixing terms of fl20l) and PUI) that the amplitudes behave differently with respect to 
the temperature, and that this provides a further method of distinguishing between the two 



sources of oscillation (see 
magnetic breakdown). 



26( 1 for the relevence of this point in the context of the theory of 



17 




c) d) 

FIG. 6: Figures showing the details of the I/2 + /1I amphtudes of {m2/mi) values of a) 0.9, b) 0.75, 
c) 0.6 and d) 0.4 . 

Our calculations can also be applied to quasi-three dimensional metals, where t is less than 
or comparable to fi. Figure [7] shows the Fourier transform of azz when fi = 5t, {001/002) = 
{1112/ nil) = ■T'o^ksT = 2.5 x Bq = rriit / {200he) , and we have set to be zero, in 

order to clarify the harmonic behavior of the amplitudes. In this case our results correspond 
to a Fourier series truncated after the second harmonics of each band, and we can see that 
the mixing effects due to scattering or chemical potential oscillations are suppressed in this 
limit. We can observe that the ratio between the first and second harmonics is ~ 2^/^ in 
the cases where we have scattering, and ~ 2^/^ in the case of oscillating /i alone - this 
difference is due to our neglecting the contribution from scattering to the oscillations in 
the latter case, which, as we can see from a comparison of the amplitudes in the graphs, 
is the dominant source of oscillations in three dimensions (as one might expect from the 

n 

results of calculations in 3D metals [27[). Naturally, at finite temperatures, the amplitudes 
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FIG. 7: Fourier transforms of the SdH oscillations in the large t limit for a window of .9 < 
Bq/B < .95 where the frequency is given in terms of B/Bq and fi = 5t, {uJi/i02) = (m2/mi) = 
.75,kBT = 2.5 X lO^^t, Bq = mit/200he, which is around 4 Tesia when mi = rUe and t = .leV 
, and we have set To to be zero. The cosine and sine functions in the amplitudes have split the 
peaks of each harmonic in two; these are labeled / and F corresponding to the 'neck' and 'belly' 
frequencies repectively. The y axis in both graphs has been rescaled by a factor of 1000; note that 
the amplitudes in the case of oscillating /i alone (A) are much smaller than those of the other 
cases. The legend in the Fourier Transform plots should be interpreted as follows: A - closed 
system, fixed r; B - open system, inter-band scattering; C - closed system, intra-band scattering; 
D - closed system, interband scattering; E - open system, intraband scattering. The plots of the 
damping factors show how the first and second order damping terms vary with their respective 
temperatures. 
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will decay correspondingly more quickly due to the effect of the temperature and Dingle 
damping factors, which can be extrapolated from their respective plots displayed in the 
figure. 

V. CONCLUSION 

We have examined possible sources of frequency mixing in the SdH oscillations of multi- 
band quasi-2D metals in the canonical and grand canonical metallic ensembles in the in- 
termediate values of the field where int > Ua- We considered a closed system with a fixed 
relaxation time r, an open system with oscillations in r arising from intra-band scattering, 
a closed system with oscillations in r arising from intra-band scattering, an open system 
with oscillations in r arising from inter-band scattering, and lastly a closed system with 
inter-band scattering with oscillations in r arising from inter-band scattering. 

In all cases other than the second we discover some measure of frequency mixing. How- 
ever, the behavior of the mixing amplitudes are slightly different for each case where they oc- 
cured, which may allow the mechanisms involved in experimental systems to be distingished. 
In the fourth and fifth cases, we observe frequency mixing due to scattering, which is inter- 
esting as we would expect such an effect in the dHvA oscillations to be virtually negligible 
gued in Section llllj) . 
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